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Abstract 

Background: Thousands of different long non-coding RNAs are expressed during embryonic development, but the 
function of these molecules remains largely unexplored. 

Results: Here we characterize the expression and function of Six30S, a long non-coding RNA that is transcribed 
from the distal promoter region of the gene encoding the homeodomain transcription factor Six3. Overexpression 
and knockdown analysis of Six30S reveals that it plays an essential role in regulating retinal cell specification. We 
further observe that Six30S regulates Six3 activity in developing retina, but does not do so by modulating Six3 
expression. Finally, we show that Six30S binds directly to Ezh2 and Eya family members, indicating that Six30S can 
act as a molecular scaffold to recruit histone modification enzymes to Six3 target genes. 

Conclusions: Our findings demonstrate a novel mechanism by which promoter-associated long non-coding RNAs 
can modulate the activity of their associated protein coding genes, and highlight the importance of this diverse 
class of molecules in the control of neural development. 



regulate mRNA transcription, acting in cis to regulate het- 
erochromatin formation at nearby genomic loci. The Xistl 
Tsix transcripts mediate X-inactivation in placental mam- 
mals [8], and Kcnqlot is important for silencing of the 
Kcnq locus resulting from parental imprinting [9]. Other 
IncRNAs regulate transcription of genes that are located 
great distances away from their own genomic loci. One 
notable example of such a trans-acting IncRNA is 
HOT AIR, which is transcribed from within specific Hox 
gene clusters, but which regulates the expression of Hox 
genes located on different chromosomes [10,11]. HOTAIR, 
Kcnqlot and Xist all mediate their effects by interacting 
with the Polycomb-repressive complex 2 (PRC2) compo- 
nent Ezh2 (enhancer of zeste homolog 2 (Drosophila)) and 
modulating histone methylation [9,11,12]. Finally, a small 
number of IncRNAs also directly interact with transcrip- 
tion factors, and potentially function as transcriptional 
coregulators [13-15]. Although the emerging picture sug- 
gests that IncRNAs may play an important and widespread 
role in regulating mammalian gene expression, a central 
and still unresolved question is how IncRNAs act in trans 
to regulate expression of specific target genes without the 
use of homologous base paring. 

O© 201 1 Rapicavoli et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BiolVlGCl Ccntrsl Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Background 

It has recently become clear that long non-coding RNAs 
(IncRNAs) comprise a large fraction of the mammalian 
transcriptome [1]. Much effort has been focused on func- 
tional analysis of IncRNAs that are processed into short 
fragments, such as microRNAs, that regulate expression 
of protein coding genes via homologous base pairing. 
However, several thousand mammalian IncRNAs have 
been identified that span multiple kilobases in length, 
and in some cases show extensive conservation at the 
nucleotide level [2-4]. 

To date, only a small number of IncRNAs have been 
functionally characterized, although this list is growing 
rapidly. Some IncRNAs act via antisense base pairing to 
block gene expression [5-7], but many show no clear 
sequence overlap with the mRNAs of protein coding 
genes. Several of these IncRNAs are known instead to 
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A complex assortment of IncRNAs is expressed in the 
developing and mature mammalian central nervous 
system, with the cellular expression patterns of nearly 1, 
000 different IncRNAs having been previously described 
[16-18]. Many show highly specific expression in specific 
brain regions and neuronal subtypes and it has been 
speculated that these IncRNAs may play a critical role 
in generating and maintaining the great cellular com- 
plexity found in the central nervous system [19,20]. 
Although a limited number of intergenic IncRNAs have 
been found to regulate neural development, their mode 
of action remains obscure [21,22]. Mechanistic insight 
into the function of one brain-expressed IncRNAs has 
come from analysis of Evf-2, a IncRNA co-transcribed 
with the homeodomain factor Dlx6. Evf-2 modulates 
transcription of Dlx6 by recruiting DLX2 and MECP2 to 
the ultraconserved ei enhancer element that is also tran- 
scribed as part of Evf-2 itself. The transcribed domain 
containing the ei sequence is essential for Evf-2 to acti- 
vate expression of Dlx6, which has raised the possibility 
that Evf-2 might regulate Dlx6 transcription at least in 
part through the formation of a RNA-DNA hybrid; this 
hybrid may in turn facilitate binding of the ei sequence 
by DLX2 and MECP2 [14,15]. 

Recent studies have also begun to address the function 
of long non-coding opposite-strand transcripts (IncOSTs), 
which are divergently co-transcribed with a broad range of 
neuronally expressed genes. Over one-third of brain- 
expressed homeodomain genes possess an associated 
IncOST, which typically spans the promoter, but not the 
transcribed region, of the protein coding gene in question 
[23,24]. Since short promoter-associated ncRNAs can reg- 
ulate expression of nearby protein coding genes [25-27], 
this has raised the possibility that these IncOSTs might 
also act in cis to selectively regulate the expression of their 
associated protein coding gene. However, although 
IncOSTs comprise a substantial fraction of all brain- 
expressed IncRNA species, their function has yet to be 
directly investigated. 

In this study, we characterize the molecular function and 
mechanism of the IncOST Six30S. Six30S is co-expressed 
with the homeodomain factor Six3, a homologue of the 
Drosophila sine oculis gene [28,29]. Like sine oculis, Six3 
plays a critical role in mammalian eye development, regu- 
lating both early eye formation and cell specification in the 
postnatal retina [30,31]. Both Six3 and Six30S are strongly 
and selectively expressed in the developing mouse retina 
and hypothalamus [23,32,33]. We use both in vivo overex- 
pression and short hairpin RNA (shRNA)-mediated 
knockdown analysis to analyze whether gain or loss of 
function of Six30S results in altered differentiation of spe- 
cific retinal cell subtypes. We also examine whether 
Six30S acts cooperatively with Six3 to regulate retinal dif- 
ferentiation, but find that Six30S does not regulate Six3 



expression levels. Finally, we provide evidence that Six30S 
can directly bind both to known transcriptional coregula- 
tors of Six3 and to histone modification enzymes, thereby 
functioning as an RNA-based transcriptional scaffold. 
These results demonstrate the mechanism by which this 
diverse class of molecules regulates cell specification dur- 
ing development. 

Results 

Genomic organization of Six3 and Six30S 

To determine which regions of Six30S to target for func- 
tional analysis, we first examined the genomic organiza- 
tion and evolutionary conservation of this IncOST using 
publicly available cDNA and genomic sequence. As pre- 
viously reported, we found that Six30S is alternatively 
spliced in both mouse and human [23], although Six30S- 
like transcripts were not identified in any other vertebrate 
species examined (data not shown). To determine which 
sequences to target for overexpression and knockdown, 
we aligned these sequences and observed that the puta- 
tive full-length mouse cDNA BC065087 contains the 
exon sequences shared by all alternative splice forms of 
Six30S (Figure 1, light grey bars). Two of these exons are 
adjacent to alternatively spliced exons of the human 
Six30S orthologue, which lie within intronic regions of 
the mouse transcript (Figure 1, dark grey bars). This 
cDNA corresponds to the most abundant isoform in neo- 
natal retina as measured by serial analysis of gene expres- 
sion (SAGE) tag abundance [32], and also matches the 
4.5-kb isoform of Six30S, previously reported to be the 
most abundant isoform expressed in embryonic brain 
[33]. We therefore selected this cDNA for analysis by 
overexpression. 

Cellular expression pattern of Six30S 

To determine if Six3 and Six30S are co-expressed in ret- 
inal progenitors, as had been previously suggested 
[23,32,33], we performed chromogenic in situ hybridiza- 
tion (ISH) on sections of embryonic day (E)16.5 and 
postnatal day (P)0.5 retina (Figure 2A-H). We observed 
that Six3 and Six30S are co-expressed in both retinal 
progenitor cells and newly post-mitotic cells of the inner 
neuroblastic layer. Additionally, Six30S is expressed in 
the developing lens (Figure 2C), while Six3 expression is 
concentrated at the interface of the inner and outer neu- 
roblastic layers (Figure 2B) at E16.5. At P0.5, Six30S 
expression appears to be more concentrated in the outer 
neuroblastic layer when compared to Six3, which is pri- 
marily found in the inner neuroblastic layer (Figure 2H). 
The differences at the interface between these layers are 
more difficult to discern at P0.5, due to the increased 
expression of Six30S in the outer neuroblastic layer (Fig- 
ure 2E-H). To determine the subcellular expression pat- 
tern of Six30S during development, we performed 
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Figure 1 Genomic structure and evolutionary conservation of Six3 and Six30S. Schematic drawing showing conservation of human and 
mouse 5ix3 and Six30S genomic structure. Light grey bars show where the predominant Six30S mouse and human forms overlap. Dark grey 
bars indicate where mouse Six30S is adjacent to conserved regions. Conservation is plotted using the PhastCons program [47]. An asterisk 
indicates the cRNA probes that were hybridized to the protein microarray in Figure 6. 



fluorescent ISH (FISH) with a probe for Six30S at E14.5, 
at which point prominent expression in progenitor cells 
of the outer neuroblastic layer is seen (Figure 21). Inter- 
estingly, though most Six30S RNA is localized in the 
cytoplasm of retinal progenitors at this stage, some 
Six30S RNA is also associated with nuclear DNA (Figure 
21, J). This subcellular distribution has also been noted 
for the IncRNA Evf-2 [14]. Using FISH in conjunction 
with immunostaining for the cytoplasmic ribosomal 



protein S6, we observed that the Six30S transcript is also 
found in the nucleus and cytoplasm when expressed in 
transfected HeLa cells (Additional file 1). 

Overexpression and knockdown of Six3 and Six30S 

Since Six30S was coexpressed with Six3 in retinal pro- 
genitors, we hypothesized that Six30S might regulate the 
expression and/or activity of Six3 in developing retina. 
To determine if this was indeed the case, we employed 
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Figure 2 Expression and subcellular localization of S/x30S and Six3 in the developing retina. (A-H) Analysis of adjacent sections using 
chromogenic ISH for Six30S and Six3 indicates substantial overlap of Six30S and Six3 expression at E16.5 (A-D) and P0.5 (E-H). (I, J) Section 
fluorescent in situ hybridization (FISH) showing Six30S expression. At El 4.5, Six30S is localized in both the nucleus (crosshairs) and the cytoplasm 



in vivo electroporation of neonatal (P0.5) mice to deter- 
mine whether overexpression and knockdown of Six30S 
phenocopied the effects of overexpression and knock- 
down of Six3. For overexpression of Six30S, we used 
full-length cDNAs corresponding to BC065087 (Figure 1) 
cloned into the pCAG plasmid [34]. Upon transfection 
into HeLa cells, robust Six30S expression was confirmed 
by FISH (Additional file 1). To visualize electroporated 
cells, all constructs were co-electroporated with plasmids 



expressing green fluorescent protein (GFP) from the 
same CAG promoter. 

We observed that overexpression of 1 [ig of Six30S 
resulted in no significant change in any major retinal 
cell type, and section immunohistochemistry revealed 
that electroporated cells showed grossly normal mor- 
phology (Figure 3B, C, white arrowheads). However, 
when dissociated cell preparations were examined, a 
significant decrease (P < 0.05) in the fraction of 
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Figure 3 Overexpression of Six30S in developing retina inhibits changes in cell fate and photoreceptor morphology observed 
following Six3 overexpression. (A) Electroporation of CAG-5/'x3 led to an increase in the fraction of GFP-positive cells expressing the 
amacrine cell marker syntaxin and a decrease in cells expressing the rod bipolar marker PKCa, CAG-5/'x305 led to a decrease in syntaxin- 
positive cells, and co-expression led to a cell composition that was indistinguishable from CAG-GFP controls. *P < 0.05. (B-E) Section 
immunohistochemistry of retinas electroporated with CAG-GFP, CAG-Six30S, CAG-Six3 or both CAG-Six30S and CAG-S/x3. White dashed lines 
define the outer third of the outer nuclear layer (OONL). Syntaxin (red) is co-immunostained with GFP (green). (C) No obvious difference is 
observed in either amacrine cell number or morphology (white arrowheads) in retinas electroporated with CAG-Six30S relative to CAG-GFP 
controls or in the fraction of the cells in the OONL. (D) An increase in amacrine cell number and in the number of cells in the OONL and a 
decrease in outer segment length is observed in the case of Six3. (E) An increase in the number of cells in the OONL is observed, but no 
difference in amacrine cell number or outer segment length is observed from controls following co-electroporation of CAG-Six3 and CAG- 
Six30S. GCL, ganglion cell layer; INL, inner nuclear layer; ONL, outer nuclear layer; OS, outer segment. (F) Laminar position of cells within 
the OONL. Electroporation of CAG-Six3 leads to a shift of rod photoreceptor cell bodies to the OONL, and this effect is not affected by co- 
electroporation of CAG-Six30S (white dashed lines in B-E). (G) Rod photoreceptors electroporated with CAG-Six3 show substantially 
shorter outer segments and this effect is reversed by co-electroporation with CAG-5/'x305. Cell type specific markers used: rhodopsin (Rho4D2), 
rod photoreceptors; glutamine synthetase (GS), Muller glia; protein kinase C alpha (PKCa), rod bipolar cells; syntaxin (Syn), amacrine cells. 



GFP-positive cells expressing amacrine cell-specific mar- 
ker syntaxin was observed relative to cells electroporated 
with the GFP control vector at P21 (Figure 3A). Since 
Six30S overexpression did not alter the fraction of GFP- 
positive cells with amacrine-like morphology and lami- 
nar position, we conclude that Six30S does not inhibit 



amacrine cell specification. The reduction in syntaxin 
expression seen following CAG-Six30S electropora- 
tion thus reflects a quantitative reduction in syntaxin 
levels or possibly reflects reduced adhesion of Six30S- 
expressing amacrine cells following dissociation for 
immunocytochemistry. 
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We next tested the effects of overexpression of pCAG- 
Six3 using in vivo electroporation. Upon electroporation, 
expression of Six3 was confirmed by section immuno- 
histochemistry. Cell fate specification in the postnatal 
retina is highly sensitive to Six3 dosage, and both gain 
and loss of function of Six3 have been reported to result 
in similar defects in bipolar cell and photoreceptor 
development [30]. We observed that electroporation of 
1 (ig of pCAG-Six3 led to a reduction in the number of 
rod bipolar cells, and detected an increase in the frac- 
tion of GFP-positive amacrine cells by P21 (P < 0.05; 
Figure 3A, D). The length of the outer segments of 
GFP-positive rod photoreceptors was also decreased 
(P < 0.05; Figure 3D, G). Finally, we observed that the 
cell bodies of the Six3 electroporated photoreceptor 
cells were primarily located in the outer third of the 
outer nuclear layer (P < 0.05; Figure 3D, white dashed 
line, and 3F), in contrast to photoreceptors electropo- 
rated with control vector, which were distributed 
throughout the outer nuclear layer. 

We then analyzed the effects of loss of function of 
Six30S and Six3. To this end, we first tested shRNAs 



for their ability to reduce expression of endogenous Six3 
and Six30S expression by in vivo electroporation of P0.5 
retina. Reduction in expression of the target gene was 
determined by analyzing expression of either Six3 pro- 
tein or Six30S RNA dissociated GFP-positive cells. For 
both Six3 and Six30S, individual shRNA constructs 
were identified that resulted in a substantial reduction 
in the average fluorescence intensity in either Six3 pro- 
tein or Six30S RNA in GFP-positive cells (Additional 
file 2). 

A significant decrease in the fraction of protein kinase C 
a (PKCa) positive rod bipolar cells was observed at P21 
following electroporation of either 1 \ig of Six3 shRNA 
construct or 1 \ig of Six30S shRNA (P < 0.05; Figure 4A), 
reminiscent of effects of retroviral overexpression of a 
dominant-negative mutant form of Six3 [30]. The fraction 
of glutamine synthetase-positive Muller glia was also 
increased (P < 0.05; Figure 4A). Section immunohisto- 
chemistry confirmed an increase in glutamine synthetase- 
positive cells with Muller glia-like morphology, at the 
expense of cells expressing the bipolar cell marker CHX10 
and demonstrating bipolar-like morphology (Figure 4B-D, 
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Figure 4 shRNA-mediated knockdown of Six3 and Six30S in the developing retina alters bipolar cell and Muller glial development. (A) 

5ix3 knockdown and Six30S knockdown both lead to an increase in Muller glia and a decrease in rod bipolar cells as measured by 
immunostaining of dissociated electroporated retinal cells. Simultaneous knockdown of both 5ix3 and Six30S resulted in a normal number of rod 
bipolar and Muller glia, but a decrease in syntaxin (Syn)-positive amacrine cells was observed. *P < 0.05. GS, glutamine synthetase. (B-E) Section 
immunohistochemistry demonstrates a decrease in GFP-positive bipolar cells following knockdown of either Six30S or Six3, but not both in 
combination (white arrowheads). Chx10 (ceh-10 homeodomain containing homolog (C. elegans)) was used as pan-bipolar cell marker. GCL, 
ganglion cell layer; INL, inner nuclear layer; ONL, outer nuclear layer; OS, outer segment. 
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white arrowheads). To further confirm that knockdown 
experiments did indeed result in a selective loss of func- 
tion of Six30S, we also inhibited Six30S function by over- 
expression of «S/#30«S-IRES-GFP. Fusion of IRES-GFP to 
IncRNAs results in a mislocalization of IncRNAs to the 
ribosome, and can produce dominant negative phenotypes 
when overexpressed [22]. When this is performed, we 
observed a reduction in rod bipolar markers and an 
increase in Muller markers, phenocopying the effects of 
Six30S knockdown (data not shown). 

Simultaneous overexpression and knockdown of Six3 and 
Six30S demonstrates non-additive effects on retinal 
differentiation 

Having observed that loss of function of Six3 and Six30S 
resulted in similar phenotypes in developing retina, we 
next investigated whether these two genes acted coopera- 
tively or independently to regulate retinal differentiation. 
We first used in vivo electroporation to overexpress 1 [ig 
of both the Six3 and Six30S constructs simultaneously in 
P0.5 retina. We observed that co-expression of Six3 and 
Six30S resulted in a normal cell fate phenotype (Figure 
3A). The length of GFP-positive rod photoreceptor outer 
segments was indistinguishable from controls (Figure 3B, 
E, G). However, the cell bodies of the electroporated 
photoreceptor cells were primarily located in the outer 
third of the outer nuclear layer (P < 0.05; Figure 3B, E, 
white dashed line, and 3F). We thus concluded that over- 
expression of Six30S was largely able to reverse the cellu- 
lar phenotypes observed following overexpression of Six3. 

As a follow up to these experiments, we simultaneously 
reduced expression of Six3 and Six30S by co-electropor- 
ating 1 [ig of the Six3 shRNA construct and 1 [ig of 
Six30S shRNA construct. We observed that simulta- 
neous knockdown of both Six3 and Six30S resulted in a 
normal number of bipolar cells and Muller glia. However, 
we observed a significant decrease in the fraction of GFP- 
positive amacrine cells, which was not observed following 
electroporation of either the Six3 or Six30S shRNA 
alone (P < 0.05; Figure 4A, E). These data show a non- 
additive effect of simultaneous loss of function of both 
Six3 and Six30S, and show that these genes functionally 
interact to regulate retinal differentiation. 

Lastly, we directly investigated the functional relation- 
ship between Six3 and Six30S in vivo by simultaneously 
overexpressing Six3 and knocking down Six30S, and vice 
versa. When we co-electroporated 1 [ig of CAG-Six3 with 
1 [ig of Six30S shRNA, we found that Six3 overexpression 
fully rescued the Six30S knockdown phenotype. Further- 
more, loss of function of Six30S fully eliminated the 
effects of Six3 overexpression. No differences in the com- 
position of GFP-positive cells were observed when com- 
pared to controls (Figure 5A). Both the morphology and 
distribution of the cells within the retina were normal 



(Figure 5C, D; Additional file 3A), and the rod photorecep- 
tor outer segment length was normal (Additional file 3B). 
We conclude that Six30S is required for many of the 
effects on retinal differentiation that result from Six3 
overexpression. 

In contrast, when 1 [ig of CAG-Six30S was co-elec- 
troporated with 1 [ig of Six3 shRNA, we observed an 
increase in glutamine synthetase-positive cells with Mul- 
ler glia-like morphology (P < 0.05; Figure 5B) and a 
decrease in PKCa-positive and CHXIO-positive bipolar 
cells (P < 0.05; Figure 5B, C, E) when compared to con- 
trols. Additionally, a decrease in cells expressing the 
amacrine cell-specific marker syntaxin was observed 
(P < 0.05; Figure 5B). These retinas, however, do not 
show any change in the number of GFP-positive cells 
with amacrine cell-like morphology. We therefore con- 
clude that, as is the case when Six30S alone was overex- 
pressed, Six30S overexpression in conjunction with 
knockdown of Six3 reduces syntaxin expression without 
otherwise affecting amacrine cell differentiation (Figure 
5B, E). The phenotype seen here is the sum of the phe- 
notypes seen following Six30S overexpression and Six3 
knockdown, and stands in contrast to the non-additive 
phenotype seen when Six30S and Six3 are both simulta- 
neously knocked down (Additional file 4). We thus con- 
clude that the phenotype observed following Six30S 
overexpression is not affected by Six3 loss of function, 
and that Six30S may also regulate syntaxin expression 
through a «S7x3-independent mechanism. 

Six30S does not directly regulate Six3 protein levels 

Although these data suggest that Six30S regulates Six3 
function in developing retina, they do not directly address 
the molecular mechanism by which Six30S is able to do 
so. A number of studies have suggested that both promo- 
ter and enhancer-associated IncRNAs can act in cis to reg- 
ulate expression levels of nearby protein-coding genes 
[14,25-27,35,36]. To investigate whether Six30S might be 
regulating Six3 function by a similar mechanism, we over- 
expressed and knocked down expression of Six30S at P0.5 
in retina, and quantified the level of Six3 protein in GFP- 
positive cells at P5.5. We did not observe any change in 
Six3 protein levels following either overexpression or 
knockdown of Six30S (Additional file 5 A, B). Further- 
more, we tested whether Six30S could regulate Six3- 
dependent regulation of a reporter construct derived from 
the sequence of the mouse Six3 promoter, which was pre- 
viously shown to be negatively regulated by Six3 [30]. To 
determine if Six30S modulates Six3-mediated auto-repres- 
sion, CAG-Six30S was cotransfected with the CMV-S/x3 
and the Six3-pro luciferase vector. As previously reported, 
overexpression of Six3 represses this reporter construct 
when tested in NIH3T3 cells. However, Six30S overex- 
pression had no effect on Six3-mediated auto-repression 



Rapicavoli et al. Neural Development 201 1, 6:32 
http://www.neuraldevelopment.eom/content/6/1/32 



Page 8 of 14 




Figure 5 Six3 overexpression rescues the Six30S knockdown phenotype whereas Six30S overexpression is not sufficient to rescue the 
loss of Six3. (A) Six30S knockdown combined with Six3 overexpression led to a cell type distribution indistinguishable from cells electroporated 
with control constructs when tested by immunostaining of dissociated electroporated retinal cells. (B) 5ix3 knockdown combined with Six30S 
overexpression led to an increase in Muller glia and a decrease in rod bipolar cells and a decrease in syntaxin (Syn)-positive amacrine cells. GS, 
glutamine synthetase. (C-E) Section immunohistochemistry of retinas confirm a distribution of GFP-positive cells that is indistinguishable from 
controls when Six30S knockdown is combined with Six3 overexpression, while Six3 knockdown and Six30S overexpression led to a decrease in 
bipolar cells (F, white arrowheads). GCL, ganglion cell layer; INL, inner nuclear layer; ONL, outer nuclear layer; OS, outer segment. 



(Additional file 6). Although 14 additional human and 
mouse cell lines were tested, no «S/x30«S-dependent effects 
on Six3-mediated autorepression could be detected (data 
not shown). We conclude that transcription of Six3 
mRNA is not regulated by Six30S. Instead Six30S RNA 
likely modulates the ability of Six3 protein to activate or 
repress expression of its target genes in developing retinal 
cells. 

Ezh2 and Eyal/3/4 directly bind Six30S 

To identify proteins that directly bind Six30S, Cy5- 
labeled Six30S RNA was used to probe human protein 
microarrays [37]. Mouse Six30S and three human splice 
variants of Six30S were transcribed in vitro and labeled 
with Cy5, and antisense cRNAs for these same tran- 
scripts were hybridized in parallel as negative controls. 
A total of five proteins that specifically and selectively 
interact with both the mouse and human forms of 
Six30S were thus identified. We found that both human 
and mouse Six30S selectively interact with Eyal, 
encoded by a homologue of the eyes absent gene of Dro- 
sophila. In addition, Six30S directly binds to the 



chromatin remodeling enzyme subunits Smarcel and 
Ezh2, as well as Enol and Ppp5c (Additional file 7). 
Members of the Eya protein family have been previously 
shown to interact with Six family proteins [38], but Eyal 
is not expressed in the mouse neuroretina [32]. There- 
fore, we also performed RNA immunoprecipitation 
(RIP) experiments on Eya3 and Eya4, which are both 
robustly expressed in embryonic retina. 

To test whether Ezh2 and Eya family members inter- 
act with mouse Six30S in vivo, co-immunoprecipitation 
experiments were performed in HEK 293T cells overex- 
pressing V5-tagged Eyal, Eya3, Eya4 and Ezh2. Expres- 
sion and immunoprecipitation of V5-tagged protein was 
verified by immunoblot (Additional file 8). We con- 
firmed that mouse Six30S RNA selectively interacts 
with each of these proteins using RIP (Figure 6), and 
conclude that these proteins interact with Six30S in 
transfected HEK 293T cells. To confirm that Six30S 
interacts with Eya family members and Ezh2 specifically, 
we examined whether Six30S interacts directly with 
Six3 protein by RIP and found that Six30S and Six3 do 
not interact in transfected HEK 293T cells (data not 
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Six30S Six30S + Six30S + Six30S + Six30S + 
EYA1 EYA3 EYA4 EZH2 



Figure 6 S/x30S-protein interactions detected via protein 
microarray. aV5-RIP specifically retrieves Eya1-V5, Eya3-V5, Eya4-V5 
and Ezh2-V5 bound to Six30S RNA but not U1 small nuclear RNA. 
Mean ± standard deviation, n = 3, relative to input is shown. 

V J 



shown). This suggests that Six30S regulates Six3 func- 
tion by facilitating interaction between Eya proteins and 
chromatin-modifying enzyme complexes. 

Discussion 

Six30S and Six3 modulate retinal cell specification and 
differentiation 

Our findings demonstrate that both gain and loss of 
function of Six30S and Six3 affect retinal differentiation 
(Additional file 4). Knockdown of Six30S resulted in a 
decrease in the fraction of bipolar cells and an increase in 
Muller glia in electroporated cells. An identical pheno- 
type was observed following knockdown of Six3 (Figure 
4). Knockdown of either Six30S or Six3 phenocopies the 
effect of retroviral overexpression of a dominant-negative 
form of Six3 in neonatal retina, where a substantial 
decrease in bipolar cells was also seen [30] . 

In contrast, overexpression of either Six3 or Six30S 
alone produced very different phenotypes. Six30S over- 
expression resulted in a reduction in syntaxin staining 
in electroporated amacrine cells, but did not alter the 
morphology or laminar position of electroporated cells. 
On the other hand, Six3 overexpression, like Six3 loss of 
function, resulted in a dramatic reduction in the number 
of bipolar cells (Figure 3A), and also resulted in 
decreased rod photoreceptor outer segment length at 
P21 (Figure 3D, G). These effects were previously 
reported following retroviral-mediated overexpression of 
Six3 [30]. Interestingly, we also observed phenotypes 
that were not reported in this previous study, such as an 
increase in amacrine cells. A number of technical differ- 
ences between these two experiments may account for 
these differences, most notably the fact that electropora- 
tion-mediated overexpression typically results in much 
higher levels of construct expression in developing 
retina than does retroviral transduction [34]. 



Six30S regulates Six3 activity in developing retina, but 
does not regulate Six3 expression levels 

If Six3 and Six30S act independently to control retinal 
cell specification, we would expect to observe purely 
additive phenotypes when their activities are both 
altered. Strikingly, when both Six30S and Six3 were 
overexpressed simultaneously, we no longer observed 
the cellular phenotypes that resulted from overexpres- 
sion of Six3 alone. Simultaneous knockdown of both 
Six30S and Six3 also showed non-additive effects. The 
loss of bipolar cells and increase in Muller glia, which 
was seen following knockdown of each gene individually, 
were not observed when expression of both Six3 and 
Six30S was reduced. Instead, a decrease in amacrine 
cells was observed, which was a phenotype not observed 
when each individual gene was knocked down. These 
non-additive phenotypes demonstrate an epistatic rela- 
tionship between Six3 and Six30S, and imply that they 
functionally interact to regulate retinal differentiation. 

Having demonstrated that Six3 and Six30S interact to 
regulate retinal differentiation, we next investigated 
whether defects in Six3 could be rescued by overexpres- 
sion of Six30S, and vice versa. When Six3 is overexpressed 
but Six30S is knocked down, none of the phenotypes nor- 
mally seen following either Six3 overexpression alone or 
Six30S knockdown alone were observed (Figure 5 A, C-D). 
These data indicate that loss of function of Six30S sup- 
presses the effects of both gain and loss of function of 
Six3. However, simultaneous overexpression of Six30S 
and knockdown of Six3 resulted in a purely additive phe- 
notype that was the sum of the effects of Six3 knockdown 
and Six30S overexpression. These findings suggest that 
Six30S acts to regulate Six3 activity, since overexpression 
of Six3 can rescue the effects of Six30S knockdown, but 
not vice versa. As Six30S overexpression leads to a reduc- 
tion in syntaxin expression even following Six3 knock- 
down, this implies that Six30S may also regulate retinal 
differentiation through a Six3-independent mechanism, 
although this remains to be characterized further. 

Since both promoter- and enhancer-associated 
ncRNAs have been previously thought to act in cis to 
selectively regulate expression of nearby protein coding 
genes [36,39], our finding that Six30S did not regulate 
Six3 expression came as something of a surprise. Our 
finding that Six30S acts in trans to regulate retinal differ- 
entiation was supported by several lines of evidence. 
First, by the finding that the activity of Six3 expressed 
from the CAG promoter can be modulated by overex- 
pression or knockdown of Six30S (Figure 5). Second, we 
have directly demonstrated that neither overexpression 
nor knockdown of Six30S has any detectable effect on 
Six3 expression in retina, nor does Six30S have an effect 
on Six3-dependent autorepression when measured by 
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luciferase analysis in transfected tissue culture cells. 
These data imply that Six30S and perhaps other IncOSTs 
do not regulate transcription of their associated protein 
coding genes, but instead act in trans to regulate the 
activity directly via protein-RNA interactions. This 
mechanism may not hold true for all promoter-associated 
IncRNAs, however, and each will have to be characterized 
independently to determine its mechanism of action. 

Six30S selectively interacts with Eya family members and 
Ezh2 

While these data demonstrate a functional relationship 
between Six30S and Six3, they still leave the precise 
mechanism by which Six30S regulates Six3 activity unre- 
solved. Our finding that both human and mouse isoforms 
of Six30S interact with multiple different members of the 
eyes absent protein family demonstrates how this might 
occur in vivo. The eyes absent {eya) gene and its mamma- 
lian homologues encode protein tyrosine phosphatases 
that function as transcriptional coregulators. Eya binds 
directly to sine oculis, and acts in conjunction with sine 
oculis in controlling eye field specification in Drosophila 
[28,29]. We also observe that Six30S directly binds the 
PRC2 subunit Ezh2, and also possibly the BAF57 subunit 
Smarcel. Additional confirmation that Six30S interacts 
with Ezh2 in vivo comes from a recent study of poly- 
comb-associated RNAs, which found that Six30S was co- 
precipitated with PRC2 in embryonic stem cells [40]. We 
propose that Six30S can modulate the expression of Six3 
target genes without necessarily regulating expression of 
Six3 itself by acting as a transcriptional scaffold. When 
Six3 directs the Six3-Eya complex to bind to specific 
genomic target sequences, Six30S may thus act as a tran- 
scriptional scaffold, recruiting histone modifying enzyme 
complexes to regulate expression of Six3 target genes. 
Although trans-acting IncRNAs have previously been 
found to regulate gene expression through recruitment of 
the Ezh2-containing PRC2 histone methyltransferase 
complex, their action is restricted to a small subset of tar- 
get genes, and the mechanism by which this occurs is 
unknown [11]. Our findings point towards a plausible 
mechanism by which this could occur. However, in the 
absence of RNA-IP data that clearly demonstrate interac- 
tion of Six30S with both Eya proteins and Six3 in devel- 
oping retina, other potential mechanisms of action of 
Six30S remain plausible. 

Several unresolved questions remain, the most obvious 
being which Eya subtype is the most relevant target of 
Six30S in vivo. Eyal is known to bind Six3 both in vitro 
and in transfected tissue culture cells [41], although the 
co-expression of Six3 and Eyal is limited to the ciliary 
margin and lens [42]. Eya2, Eya3 and Eya4, however, are 
all co-expressed with Six3 and Six30S in both the devel- 
oping retina and in other forebrain regions [42]. Eya4 



and Six3 have been demonstrated to interact in a study 
of holoprosencephaly, indicating that Eya4 may be the 
functional bridge between Six3 and Six30S in the devel- 
oping ventral forebrain [43]. Our findings indicate that 
Six30S may also regulate retinal differentiation by Six3- 
independent mechanisms. Since Eya family members can 
interact with other classes of transcription factors, 
including Hox and Tlx family members [44], this sug- 
gests that Six30S might also modulate the expression of 
genes regulated by these transcription factors. 

Many homeodomain transcription factors that are 
essential for central nervous system development possess 
an associated IncOST, including Pax6, Rax, Vax2 and Six6 
[23,45]. Each of these proteins exhibits complex and 
dynamic expression during development, and their expres- 
sion patterns often diverge considerably from those of 
their associated IncOST. Although their target sequences 
are present in all cells, the genes that are directly regulated 
by each protein can vary considerably, both at different 
developmental stages and in different cell types. Models 
that explain how this context-dependent regulation of 
transcription factor activity occurs have typically empha- 
sized combinatorial regulation by other proteins, or cell- 
and stage-specific epigenetic modifications. Our finding 
that IncOSTs can modulate the activity of their associated 
transcription factors lends an additional layer of complex- 
ity to these models, and suggests that this diverse class of 
molecules may play a critical function in generation of cell 
subtype diversity within the developing central nervous 
system. 

Materials and methods 

Animals 

Pregnant CD-I mice were purchased from Charles River 
Breeding Laboratories, Wilmington, MA, USA. Animal 
experiments were approved by the Johns Hopkins Uni- 
versity Institutional Animal Care and Use Committee. 

In vivo electroporation 

In vivo electroporation of mouse retina was performed as 
described [34]. Retinas were electroporated at P0.5. Dis- 
sociated cells and section immunohistochemistry data 
shown here were performed at P21. Error bars represent 
± standard error of the mean for at least three indepen- 
dent electroporated retinas. A two-tailed Students £-test 
was performed to determine the P-value. Control and 
experimental constructs were co-electroporated with 0.2 
\ig of pCAG-GFP to readily visualize electroporated cells 
by GFP [34]. For overexpression experiments, 1 \ig of 
pCAG-GFP was injected with either 1 [ig of pCAG-Six3 
or 1 \ig of CAG-Six30S. For simultaneous overexpression 
of Six3 and Six30S, retinas were electroporated with 1 \ig 
of pCAG-Six3 and 1 \ig of CAG-Six30S. For knockdown 
of Six3 and Six30S, 1 \ig of U6-control was injected with 
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either 1 [ig of U6-Six3 or 1 [ig of U6-Six30S. For simulta- 
neous knockdown, 1 [ig of U6-Six3 and 1 [ig of U6- 
Six30S were injected. For overexpression of Six3 in com- 
bination with knockdown of Six30S, retinas were electro- 
porated with 1 \ig of pCAG-Six3 and 1 \ig of U6-Six30S 
were injected. Finally, for overexpression of Six30S 
in combination with knockdown of Six3> 1 [ig of pCAG- 
Six30S and 1 [ig of U6-«S7;\;3 were injected and 
electroporated. 

In situ hybridization 

For cryosections, untimed E14.5 and E16.5 embryos and 
P0.5 neonates were dissected, fixed overnight in 4% par- 
aformaldehyde in PBS at 4°C and cryoprotected overnight 
in 30% sucrose/PBS at 4°C before being embedded in 
OCT compound (Sakura, Torrance, CA, USA) on dry ice. 
Cryosections (20 [im) were cut on a cryostat. Section ISH 
methodology was as previously described [32] with the 
exception that probe BC065087 was used and a Tyramide 
Signal Amplification system (TSA Plus, Perkin Elmer, 
NEL 744, Waltham, MA, USA) combined with an antidi- 
goxigenin-HRP antibody (1:1, 000; Roche, Indianapolis, 
IN, USA ), for fluorescent ISH. Fluorescent ISH sections 
were counterstained with DAPI. Fluorescent ISH samples 
were photographed on a Zeiss Meta 510 confocal micro- 
scope. Chromogenic images were visualized on a Zeiss 
Axioskop2 microscope. 

Immunohistochemistry 

For cryosections, electroporated eyes were harvested 21 
days after electroporation, fixed for 1 hour in 4% parafor- 
maldehyde in PBS at 4°C and cryoprotected overnight in 
30% sucrose/PBS at 4°C before being embedded in OCT 
compound (Sakura) on dry ice. Cryosections (20 (im) were 
cut on a cryostat. Retinal cryosections were immunos- 
tained as described except that sections were post fixed in 
4% paraformaldehyde for 5 minutes prior to blocking [21]. 
Samples were photographed on a Zeiss Meta 510 confocal 
microscope. 

For cell marker analysis, dissociation and immunos- 
taining were performed as described [34] except that reti- 
nas were harvested at P21. Samples were visualized and 
quantified on a Zeiss Axioskop2 microscope. To quantify 
the effects of Six3 knockdown, retinas were electropo- 
rated with Six3 shRNA, harvested at P4.5, dissociated, 
and immunostained as described above with anti-Six3 
and detection with Alexa568 goat anti-guinea pig IgG 
and rabbit anti-GFP and detection with Alexa488 goat 
anti-rabbit. For Six3 quantification, Six30S overexpres- 
sion and SixOS shRNA electroporated retinas were har- 
vested at P5.5, dissociated and immunostained as 
described above with anti-Six3 and detection with 
Alexa568 goat anti-guinea pig IgG and rabbit anti-GFP 
and detection with Alexa488 goat anti-rabbit. Samples 



were visualized on a Zeiss Axioskop2 microscope and 
signal intensity was quantified with Velocity 4.0 (Perkin 
Elmer) software by calculating the average signal inten- 
sity per cell and normalized to cell size. 

In situ hybridization and immunocytochemistry 

For Six30S knockdown quantification, retinas were har- 
vested at P4.5 and dissociated as described above. FISH 
was performed as described using TSA Plus Cyanine3 
kit (1:125; Perkin Elmer, NEL 744) followed by staining 
as described (with rabbit anti-GFP and detection with 
Alexa488 goat anti rabbit). Samples were visualized and 
quantified as above. For all dissociations, nuclear DNA 
was visualized with DAPI counterstaining. Cell counts 
were analyzed using the two-tailed Student's £-test. A 
minimum of three retinas were counted for each con- 
struct examined using dissociated immunocytochemis- 
try, with 100 to 300 GFP-positive cells per retina 
counted for each marker tested. 

HeLa cell FISH followed by immunocytochemistry was 
performed as previously described [22] except that 20 
cells were counted. 

Luciferase reporter assays 

NIH 3T3 cells were grown in DMEM supplemented with 
10% FCS. Cells were transfected with Fugene6 (Roche) per 
the manufacturer's instructions. Cells were transfected 
with 500 ng of luciferase reporter construct and 50 ng of 
the expression plasmids for Six3 and Six30S. Cells were 
harvested 2 days post-transfection. Luciferase was mea- 
sured per manufacturers instructions with the Dual-Luci- 
ferase Reporter System (Promega, E1910, Madison, WI, 
USA). pTK-Renilla (50 ng) was co-transfected to control 
for transfection efficiency. 

Synthesis of RNA probes for protein microarray 

BC065087 (5 [ig) was linearized with Dral, and 5 [ig of 
BM663835, BX115070, and BM451513 were linearized 
with Notl. Linearized DNA was then phenol-chloroform 
extracted, ethanol precipitated and resuspended in 10 [i\ 
water. Probes were in vitro transcribed with T7 poly- 
merase per manufacturer's instructions using the Ribop- 
robe Combination system kit (Promega, P1405), spiking 
the reaction with 1 [i\ of 5 mM Cy5 labeled CTP (GE 
Healthcare, 25-8010-87, Piscataway, NJ, USA). The RNA 
probe was ethanol precipitated with LiCl and resus- 
pended in 1 mM EDTA. 

Protein microarray analysis 

Transcription factor/RNA binding protein chips were 
generated as previously described [37]. Protein chips 
were preblocked with Superblock Buffer (Thermo 
Scientific, 37516, Waltham, MA, USA) supplemented 
with 10 [ig ml" 1 sperm DNA, 2 mM MgCl 2 , 2 mg ml" 1 
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BSA for 1 hour at room temperature. RNA was hybri- 
dized at 250 nM concentration in binding buffer (PBS 
supplemented with 2 mM MgCl 2 , 2 mg ml" 1 BSA) at 
room temperature for 1 hour. The slides were washed 
four times in TBST, dried and scanned by a GenePix 
400B scanner. Data were analyzed with GenePix Pro 
6.10 as previously described [46]. 

RNA immunoprecipitation experiments 

HEK 293T cells were grown in DMEM supplemented 
with 10% FCS. Cells were transfected with Fugene6 
(Roche) per the manufacturer's instructions. Ten million 
cells were transfected 24 hours post-plating with 5 [ig of 
pCAGIG-V5, pCAGIG-Six3-V5, pCAGIG Eyal-V5, 
pCAGIG Eya3-V5, pCAGIG Eya4-V5 or pCAGIG Ezh2- 
V5, and pCAG-Six30S. Cells were harvested 2 days 
post-transfection, lysed and precipitated essentially as 
previously described [11], except that 5 [ig anti-V5 anti- 
body was used (1:5, 000; Invitrogen, R96025, Carlsbad, 
CA, USA) with 50 \i\ protein A-agarose (Invitrogen, 
15918-014) in each immunoprecipitation reaction. For 
each sample, 10% of total volume was set aside after 
lysis for RNA extraction and 5% set aside for immuno- 
blot analysis. After precipitation, 80% of the beads were 
resuspended in Trizol and 20% were resuspended in 
Laemmle buffer. RNA was Trizol extracted and resus- 
pended in 50 \i\ of nuclease free water. RNA was then 
DNAse treated with Turbo DNA-free (Ambion, 
AM1907, Carlsbad, CA, USA) per manufacturer's 
instructions. RNA was quantified by quantitative RT- 
PCR using Brilliant II SYBR Green QRT_PCR Master 
Mix (Agilent, 600825, Santa Clara, CA, USA) on a 
Roche LightCycler480. No-RT controls were simulta- 
neously performed to demonstrate that signal was not 
from DNA contamination. 

Expression of V5-tagged protein was confirmed by 
western analysis using anti-V5 antibody (1:10, 000; Invi- 
trogen, R96025) dilution in 5% milk in PBST, and 
detected with horse radish peroxidase goat anti-mouse 
IgG (1:10, 000; Santa Cruz, sc-2031, Santa Cruz, CA, 
USA) and ECL Western Blotting Detection System (GE 
Healthcare, RPN2132 Piscataway, NJ, USA) per the man- 
ufacturers' instructions. 

Full details of all plasmids and antibodies used in this 
study are included in Additional File 9. 

Additional material 



Additional file 1: Six30S is localized equally in the nucleus and 
cytoplasm. HeLa cells were transfected with Six30S constructs and RNA 
location was analyzed by FISH followed by immunohistochemistry 
against the cytoplasmic S6 ribosomal protein. Cytoplasmic Six30S RNA 
was identified by localization with S6 protein. The relative proportion of 
nuclear Six30S, defined as FISH signal that did not colocalize with S6 
protein, is indicated. N = 20 cells. 



Additional file 2: Confirmation of shRNA-mediated knockdown of 
endogenous Six3 and Six30S in developing retina. (A E) A construct 
encoding either control hairpin, 5ix3 targeted hairpin or Six30S targeted 
hairpin was electroporated in vivo, into P0.5 retina and harvested at P4.5, 
and dissociated, and immunocytochemistry against Six3 and GFP, or FISH 
against Six30S was performed followed by immunostaining for GFP. GFP- 
positive cells were counted to analyze the fluorescence intensity for each 
cell that expressed (A) Six3 or (D) Six30S. At least 100 GFP-positive cells 
from three different electroporated retinas were counted for each 
construct tested. Error bars represent standard error for at least three 
independent retinas. (A) P < 0.05; (D) P < 0.05. (B-E) Examples of 
dissociated cells positive for GFP and either Six3 protein or Six30S RNA. 

Additional file 3: Six30S knockdown rescues changes in retinal cell 
morphology and position observed following overexpression of 

Six3. (A, B) Six30S knockdown combined with Six3 overexpression 
rescues the effects observed on rod photoreceptor cell body position (A) 
and on photoreceptor outer segment length observed with 5ix3 
overexpression (B). 

Additional file 4: Summary of overexpression and knockdown data 

for Six3 and Six30S. These results demonstrate that co-expression of 
Six30S and 5ix3 rescues the phenotypes observed with Six3 
overexpression except that the photoreceptors are displaced in the outer 
third of the outer nuclear layer. Simultaneous knockdown of Six3 OS and 
Six3 results in a novel phenotype, fewer amacrine cells. Additionally, 
Six30S overexpression rescues the phenotype of knockdown of 5ix3. 
However, expression of Six3 combined with knockdown of Six30S results 
in an additive phenotype. 

Additional file 5: Cellular levels of Six3 protein are not affected by 
overexpression or knockdown of Six30S. (A) Overexpression of Six30S 
at P0.5 does not affect Six3 protein levels when measured at P5.5. (B) 
shRNA-mediated knockdown at P0.5 of Six30S does not affect Six3 
protein levels when measured at P5.5 

Additional file 6: Autorepression of Six3 is not affected by Six30S. 

Transfection of CMV-S/'x3 together with the Six3-pro luciferase reporter 
into NIH 3T3 led to > 70% reduction in luciferase expression. Co- 
transfection of CMV-S/x3 and CAG-Six30S with the 5/x3-pro luciferase 
reporter was not significantly different than CMV-5/x3 alone. 

Additional file 7: List of proteins that interact with the mouse and 
human forms of Six30S from the transcription factor/RNA binding 
protein microarray. 

Additional file 8: V5 expression constructs are expressed in HEK 
293T cells. aV5 western analysis demonstrating that Eya1-V5, Eya3-V5, 
Eya4-V5 and Ezh2-V5 are expressed. 

Additional file 9: Supporting Information. Includes plasmid 
information, shRNA sequences and antibodies used. 
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BSA: bovine serum albumin; DMEM: Dulbecco's modified Eagle's medium; E: 
embryonic day; FCS: fetal calf serum; FISH: fluorescent in situ hybridization; 
GFP: green fluorescent protein; ISH: in situ hybridization; IncOST: long non- 
coding opposite-strand transcript; IncRNA: long non-coding RNA; ncRNA: 
non-coding RNA; P: postnatal day; PBS: phosphate-buffered saline; PKC: 
protein kinase C; PRC2: Polycomb-repressive complex 2; RIP: RNA 
immunoprecipitation; shRNA: short hairpin RNA; TBST: Tris- buffered saline 
with Tween-20. 
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